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ABSTRACT 

We present optical integral field spectroscopy of the circum-nuclear gas of the Seyfert 2 galaxy NGC 
1386. The data cover the central 7" x 9" (530 x 680 pc) at a spatial resolution of 0'.'9 (68 pc), and the 
spectral range 5700-7000 A at a resolution of 66 km s“^. The line emission is dominated by a bright 
central component, with two lobes extending 3" north and south of the nucleus. We identify three 
main kinematic components. The first has low velocity dispersion (d ~ 90 km s“^), extends over 
the whole field-of-view, and has a velocity field consistent with gas rotating in the galaxy disk. We 
interpret the lobes as resulting from photoionization of disk gas in regions where the AGN radiation 
cones intercept the disk. The second has higher velocity dispersion (d ~ 200 km s“^) and is observed 
in the inner 150 pc around the continuum peak. This component is double peaked, with redshifted 
and blueshifted components separated by ^ 500 km s“^. Together with previous HST imaging, these 
features suggest the presence of a bipolar outflow for which we estimate a mass outflow rate of M > 
0.1 Mq yr“^. The third component is revealed by velocity residuals associated with enhanced velocity 
dispersion and suggests that outflow and/or rotation is occurring approximately in the equatorial 
plane of the torus. A second system of velocity residuals may indicate the presence of streaming 
motions along dusty spirals in the disk. 

Subject headings: black hole physics — galaxies: individual: NGG 1386 — galaxies: active — galaxies: 
nuclei — galaxies: kinematics — galaxies: Seyfert. 


1. INTRODUCTION 

Gas accretion onto a supermassive black hole (SBH) is 
believed to be the mechanism which powers active galac¬ 
tic nuclei (AGNs, e.g. Lynden-Bell 1969; Begelman et al. 
1984). However, understanding how mass is transferred 
from galactic scales (kiloparsecs) down to nuclear scales 
(parsec and sub-parsecs) to create a reservoir of gas to be 
accreted has been a long-standing problem in the study 
of nuclear activity in galaxies. 

Theoretical studies and simulations have shown that 
non-axisymmetric potentials efficiently promote gas in¬ 
flow towards the inner regions of galaxies: close encoun¬ 
ters and galactic mergers have been identified as mech¬ 
anisms capable of driving gas from large scales (10 kpc) 
down to the inner kiloparsec (e.g. Hernquist 1989; Di 
Matteo et al. 2005; Springel et al. 2005b, a; Hopkins et al. 
2006). Stellar bars (and bars-within-bars) are also im¬ 
portant non-axysimmetric features deemed able to fun¬ 
nel gas from the inner kpc down to few hundreds of pc 
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from the SBH (e.g. Shlosman et al. 1989; Friedli & Mar¬ 
tinet 1993; Englmaier & Shlosman 2004). Simulations 
presented in Hopkins & Quataert (2010) suggest that, 
in gas-rich systems, on scales of ^ 10 - 100 pc, inflows 
are achieved through a system of nested gravitational in¬ 
stabilities, similar to the bars-within-bars scenario, but 
with a range of morphologies, such as nuclear spirals (in¬ 
cluding one- and three-armed spirals), bars, rings, barred 
rings, clumpy disks and streams. 

Observations support the hypothesis that large-scale 
bars funnel gas toward the center of galaxies (e.g. Gren- 
shaw et al. 2003a). Imaging has revealed that struc¬ 
tures such as small-scale disks or nuclear bars and as¬ 
sociated spiral arms are frequently observed in the inner 
kiloparsec of active galaxies (e.g Erwin & Sparke 1999; 
Pogge & Martini 2002; Laine et al. 2003). The most 
common nuclear structures are dusty spirals, estimated 
to reside in more than half of active galaxies (Martini 
et al. 2003). Strong correlation between the presence 
of nuclear dust structures (filaments, spirals and disks) 
and activity in galaxies has also been reported (Simoes 
Lopes et al. 2007). This correlation between dust struc¬ 
tures and activity, along with the enhanced frequency of 
dusty spirals, supports the hypothesis that nuclear spi¬ 
rals are a mechanism for fueling the SBH, transporting 
gas from kiloparsec scales down to within a few tens of 
parsecs of the active nucleus. 

In recent years integral field spectroscopy of nearby 
AGNs has provided a clearer view of the processes of 
AGN fueling on scales in the range ^10 - 100 pc, re¬ 
vealing kinematical features consistent with gas inflow¬ 
ing along bars and nuclear spirals (e.g. Knapen et al. 
2000; Eathi et al. 2006; Storchi-Bergmann et al. 2007; 
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Riffel et al. 2008; Muller Sanchez et al. 2009; Riffel & 
Storchi-Bergmann 2011; Schnorr Muller et al. 2011; Rif¬ 
fel et al. 2013; Combes et al. 2014; Davies et al. 2014; 
Schnorr-Miiller et al. 2014a, b). 

Removal of angular momentum can be achieved not 
only through gravitational torques, but also via the out¬ 
flows, or winds, originating from the interaction between 
ionized gas and the magnetic field (e.g. Rees 1987). Out¬ 
flows are prevalent among powerful AGNs (e.g. Cren¬ 
shaw et al. 2003b; Barbosa et al. 2009), however stud¬ 
ies of nearby Seyferts suggest that compact outflows 
100 pc in extent) with velocities of ^ 100 km s“^ 
and mass outflow rates of a few solar masses per year 
are common even in low-luminosity AGNs (e.g. Storchi- 
Bergmann et al. 2007; Davies et al. 2009; Schnorr Muller 
et al. 2011; Miiller-Sanchez et al. 2011; Davies et al. 
2014). 

Outflows transfer mechanical energy to the host inter¬ 
stellar medium, therefore quantifying mass (and kinetic 
energy) outflow rates is fundamental to understanding 
the process and pace of galaxy evolution (e.g. Di Mat- 
teo et al. 2005). Key questions are: how are the inflow 
and outflow rates related? Does the presence of one pre¬ 
clude the other? In order to answer such questions we 
have assembled a hard X-ray selected sample of 20 nearby 
(z <0.01) Seyferts covering a range of accretion power. 
With this sample we aim to (1) determine if there is a 
transition in the behavior of the circumnuclear gas kine¬ 
matics from low-luminosity to high-luminosity AGNs and 
investigate the effect of other factors such as host galaxy 
properties; (2) measure the mass inflow/outflow rate and 
relate it to AGN accretion power. 

Here we present a quantitative analysis of the gas and 
stellar kinematics in the inner 250 parsec of NGG 1386 
as derived from integral field spectroscopic (IFS) obser¬ 
vations carried out with the GEMINI South telescope. 

NGG 1386 is a Sb/c spiral galaxy (Malkan et al. 1998), 
located in the Fornax cluster (Sandage 1978), which 
hosts a Seyfert 2 nucleus. At a distance of 15.6 Mpc 
(determined from surface brightness fluctuations, Jensen 
et al. 2003) the linear scale is 76 pc arcsec“^. The in¬ 
clination with respect to the line of sight is 0 ^ 65°, as 
derived from the ratio of the apparent axes using values 
from the 2MASS catalog (Jarrett et al. 2003) available 
in NED^ 

Being one of the nearest Seyferts, NGG 1386 has been 
extensively studied. Sandage (1978) classified NGG 1386 
as a high excitation galaxy on the basis of its high 
[OIII]A5007/H/d line intensity ratio. High resolution ob¬ 
servations at 8.4 GHz show the presence of a non-variable 
compact source confined within the central arcsecond 
and elongated southward in position angle PA ^ 170 - 
175° (Nagar et al. 1999; Mundell et al. 2009). HST imag¬ 
ing of the [OHI] and Ha + [Nil] emission reveals a jet-like 
distribution of ionized gas organized in knots along the 
north-south axis and extending over a distance of about 
2" from the nucleus (Ferruit et al. 2000). Despite the 
close alignment between ionized gas and radio emission, 
direct comparison suggests that there is no detailed as- 

® NED is the NASA/IPAC Extragalactic Database (NED) which 
is operated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronautics and 
Space Administration. 


sociation amongst the two (Mundell et al. 2009). Early 
2D spectroscopic observations (Weaver et al. 1991; Rossa 
et al. 2000), supported by more recent studies (e.g. Schulz 
& Henkel 2003; Bennert et al. 2006), have already re¬ 
vealed the presence of a complex nuclear spectrum show¬ 
ing the presence of a strong narrow component, following 
rotational motion, superposed on a broader component 
suggestive of high-velocity outflows or inflows. 

The paper is organized as follows: observations and 
data reduction are presented in §2, the data analysis is 
described in §3. Our results are presented in §4 and 
discussed in §5. In §6 we propose an interpretation of the 
results summarizing our main findings and the proposed 
model in §7. 

2. OBSERVATIONS 

NGG 1386 was observed on September 23, 2011 with 
the Integral Field Unit (IFU) of the Gemini Multi-Object 
Spectrograph (GMOS; Allington-Smith et al. 2002; Hook 
et al. 2004) mounted on the Gemini South Telescope 
(program ID: GS-2011B-Q23). Two adjacent fields were 
observed, each covering 7x5 arcsec^, resulting in a total 
angular coverage of 7.5 x 9.5 arcsec^ with O'!2 sampling, 
centered on the nucleus and extended along the galaxy 
semi-major axis (position angle PA = 25°). In the top 
panels of Fig. 1 we show the observed field-of-view (FOV) 
on an HST image, the corresponding “structure map” 
created via the contrast-enhancement technique defined 
by Pogge & Martini (2002), and a continuum image ob¬ 
tained from our observation. 

In order to account for dead fibers in the IFU, ± 0?3 
spatial dithering was applied along both axes while a 5 
nm spectral dithering was applied to account for CCD 
gaps. This resulted in 12 exposures of 375 seconds each. 
From measurements on the standard star, observed im¬ 
mediately after the galaxy, we determined that the seeing 
limited the spatial resolution to 0?9 ± 0?1, corresponding 
to 68 ± 8 pc at the galaxy distance. 

In order to cover the wavelength range 5600-7000 
A, which includes the emission lines [FeVH]A6087, 
[OIJA6302, Ha+[NH]AA6548,6583, [SH]AA6716,6731 and 
several stellar absorption features, we used the IFU in 
two-slits mode with the grating GMOS R400 in combi¬ 
nation with the r(630 nm) filter, this yields a spectral 
resolution R = 1918 corresponding to a velocity disper¬ 
sion cr 66 km s“^. A portion of the spectrum including 
the most prominent emission lines is shown in the bot¬ 
tom panel of Fig.l for some representative regions of the 
FOV. 

Data reduction was performed using the IRAF^ pack¬ 
ages provided by the GEMINI Observatory and specifi¬ 
cally developed for the GMOS instruments^. The data 
reduction process includes bias, overscan and sky sub¬ 
traction, flat-fielding, trimming, wavelength and flux cal- 

^ IRAF is the Image Reduction and Analysis Facility, a general 
purpose software system for the reduction and analysis of astro¬ 
nomical data. IRAF is written and supported by the National 
Optical Astronomy Observatories (NOAO) in Tucson, Arizona. 
NO AO is operated by the Association of Universities for Research 
in Astronomy (AURA), Inc. under cooperative agreement with the 
National Science Foundation. 

http://www.gemini.edu/sciops/data-and-results/ 
processing-software?q=node/11822 
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Figure 1. Top left: HST image of NGC 1386 (WFPC2/F606W-PC, proposal ID: 5446). The vertical dashed line marks the semi¬ 
major axis (PA = 25° east of north). Top center: zoom-in on the box in the left panel showing the “structure map” (Pogge & Martini 
2002) corresponding to the observed field-of-view. The map was derived from the image in the left panel. The dashed line indicates the 
approximate position angle (PA = 76° east of north) and extension of a faint bar-like [OIII] emission observed by Permit et al. (2000). 
Top right: continuum map as obtained from our data; labels mark spaxels from which the representative spectra plotted below have been 
extracted. Bottom: representative spectra from individual spaxels A, C and N, as marked in the top right panel. 


ibration, building of the data cubes at a sampling of 0.1 
X 0.1 arcsec/pixel, final alignment and combination of 
the 12 cubes. The resulting cube contains 7030 spectra. 

The typical accuracy in the wavelength calibration is 
0.14 A (or 6 km s ^ at A = 6583 A). Flux calibration is 
expected to be accurate to ~ 10%. 


We used customized IDL^^ routines to model the 
continuum and the profiles of the most prominent 
emission lines (Ha, [Nil]AA6548,6583, [SII]AA6716,6731, 
[OIJA6302 and [FeVII]A6087). Gaussian profiles were fit¬ 
ted to the lines in order to derive centroid velocities, ve¬ 
locity dispersions and fluxes. Four separate runs were 
performed to fit the Ha + [Nil] lines, the [SH] doublet, 
the [01] line and [FeVH]. Given the limited range in wave- 


3. EMISSION LINE FITTING 


IDL, or Interactive Data Language, is a programming lan¬ 
guage used for data analysis and visualization. 
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Figure 2. Gaussian fits to the representative spectra shown in 
Fig.l. Top: three component fit. Middle: two component fit. 
Bottom: single component fit. Crosses show the data points and 
the grey solid lines show the summed fit, with dashed and dashed- 
dotted lines showing the individual components as labelled. The 
lower sub-panel in each figure shows the residuals of the fit. 


length of the fitted spectral regions, the local continuum 
adjacent to the lines was modeled with a first order poly¬ 
nomial. 

Visual inspection of the spectra suggests that, besides 
[FeVII], emission line profiles within ^ 1" of the nu¬ 
cleus mostly result from the superposition of one rel¬ 
atively narrow component (d ~ 90 km s“^) and two 
broader components (d ~ 200 km s“^), see for exam¬ 
ple the [Nil] and [SII] profile for spectrum N in Fig.l). 
Hereafter we will refer to the components as “narrow” 
and “broad” (note, however, that the broad components 
are not from the “classical” broad line region of AGN 
where line widths exceed 1000 km s“^). Beyond the cen¬ 
tral 2" the lines are clearly unblended and single-peaked 
requiring only one Gaussian to be satisfactorily modeled. 
Between 1" and 2" the line profiles mostly present a nar¬ 
row core and a broad base. We interpret the latter as 
due to the blending of the two broad components and 
we fit these lines with two Gaussians: one for the narrow 
core and one for the broad base. Note, however, that for 
a fraction of the spaxels in this region the spectra can 
still be successfully modeled with three Gaussians. 

In order to fit the HQf+[NII] lines we assumed, for each 
component of the fit (i.e. the narrow and the two broad 
components), that Ra and [NIIJA6548 have the same red- 
shift and width as [NIIJA6583 (the strongest line). There¬ 
fore, the velocity and velocity dispersion of each compo¬ 
nent of [NIIJA6548 and Ha are the same as those of the 
corresponding components of [NIIJA6583. The amplitude 
of each [NIIJA6548 component was fixed at the theoretical 
value of 1/2.96 of that of the corresponding [NIIJA6583 
amplitude. A similar approach was used to fit the [SII] 
doublet. 

A simpler double peaked profile gives a better represen¬ 
tation of the [FeVII] line which, as it lacks a detectable 
narrow component, was fitted with two Gaussians of the 
same width. Because of its high ionization potential 
{E = 125 eV, Gorliss & Sugar 1982), this line is believed 
to be produced by a combination of hard UV continuum 
and hot collisionally ionized plasma (Contini et al. 1998). 
Given the high energies required for its production, the 
[Fe VII] line is believed to originate close to the AGN (a 
few tens of pc, Ferguson et al. 1997). 

Three Gaussians fit: in order to determine the region 
over which the three Gaussian fit is necessary, we require 
that for each component the line amplitude is at least 
twice the standard deviation in the continuum adjacent 
to the line. We also assume that the two broad compo¬ 
nents have the same width (without this condition, the 
fit would be be under constrained in many spaxels, where 
the profiles are heavily blended). With this approach we 
fit three Gaussian components to each of Ha, [Nil] and 
[SII] mostly within a radius of 1" from the nucleus, but 
also at larger distances, between 1" and 2" (the nucleus, 
denoted N, is identified with the continuum brightness 
peak). Due to the lower S/N, three Gaussians were fit¬ 
ted to [01] only for a handful of pixels within ^ 077 from 
the center. An example of the three components fit is 
shown in the top panel of Fig. 2 for [Nil] and Ha. 

Two Gaussians fit: broad components are still visible 
as a broad base mainly between 1" - 2" from the nu¬ 
cleus. However, for many pixels the two components are 
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Figure 3. Map showing the regions where one, two or three Gaus- 
sians were fitted to the emission lines. No significant Hq- emission 
was detected in the darker region labelled as ■ 

blended and the S/N is insufficient to justify fitting three 
components. Therefore in this region we mostly fit two, 
representing the broad base and the narrow core, respec¬ 
tively. An example of a two-components fit is shown in 
the middle panel of Fig. 2. This fit was performed with 
the goal of extracting the narrow component. 

Single Gaussian fit: a single Gaussian fit was used be¬ 
yond the central 2" where the broad components become 
negligible and the two-components fit failed. An example 
of single component fit is shown in the bottom panel of 
Fig. 2. A map showing the number of components fitted 
to the spectra in a given spaxel is shown in Fig. 3. 

To determine accurate line intensity ratios representa¬ 
tive of the various kinematic components, five subregions 
of the FOV were identified (see velocity map for the nar¬ 
row and broad blue component in Fig. 6). The subre¬ 
gions represent, respectively, the northern and southern 
“lobes” (subregions A and B), regions associated with ve¬ 
locity residuals approximately perpendicular to the axis 
of the lobes (see §5.2 for details), east and west of the 
nucleus (G and D), and the nucleus itself (N). The spec¬ 
tra obtained for sub-regions B, D and N are presented 
in Fig. 4; the spectra for subregions A and G are very 
similar to those of subregion B and D respectively. 

The continuum and all the most prominent emission 
lines ([01], Ha, [Nil], [SII]) were fitted simultaneously 
with the constraint that all the lines are characterized by 
the same redshift and the same width as the [NII]A6583 
line. With these assumptions we fitted single Gaussians 
to the emission lines in the subregion B spectrum and two 
Gaussians to each line in subregion A. This was done to 
take into account the non-Gaussian base of the lines. In 
this case, the total flux of each line was considered to be 
the sum of the fluxes of the two components. Line emis¬ 
sion is weak relative to the stellar continuum in subre¬ 
gions C and D, therefore stellar templates were fitted and 
subtracted. The templates, a subset of the models pre¬ 
sented in Vazdekis (1999), represent an old stellar popu¬ 
lation of spectral type 0-M and metallicity -0.7 < [Fe/H] 
-kO.2. After this procedure, line profiles in subregion G 


and D were fitted with the same approach described for 
subregion A. The profiles of subregion N, which includes 
the nucleus, were fitted with three components: a nar¬ 
row core and two broader Gaussians characterized by the 
same width. In order to minimize contamination from 
Ha, the blueshifted broad component was fitted to the 
blue side of the [NH]A6548. 

The stellar continuum was not removed prior to fitting 
the subregion A, B and N spectra as this did not produce 
any significant difference in the results. 

3.1. Uncertainties 

Errors on the measured quantities were derived from 
Monte Carlo simulations. For each spaxel, we con¬ 
structed one hundred realizations of the spectrum by 
adding Gaussian noise with amplitude comparable to the 
noise measured in the original spectrum. Mean values 
and standard deviations for the centroid velocities, ve¬ 
locity dispersions and fluxes were derived for each spaxel, 
with the standard deviation of the distribution in each 
parameter being taken as the uncertainty. 

4. RESULTS 

4.1. Gas kinematics 

4.1.1. Velocity and velocity dispersion 

As [NH]A6583 has the highest S/N, we use results ob¬ 
tained from the fits to this line to represent and model 
the gas velocity field. The velocities recovered from other 
emission lines are consistent within the errors. 

We first show maps obtained by fitting a single Gaus¬ 
sian to the [Nil] line throughout the entire FOV. Cen¬ 
troid velocity and velocity dispersion maps obtained with 
this approach are presented in Fig. 5 after the subtraction 
of a systemic velocity Vsys = 786 km s“^ determined from 
the modeling of the gas kinematics (see §5.2 for details). 
The velocity map shows the characteristic pattern of a 
rotating disk upon which other components are super¬ 
posed. There are two clear deviations from the rotation 
pattern: (i) a redshifted blob located at the nucleus and 
(ii) a blueshifted region to the north-west of the nucleus. 
Multicomponent fits and visual inspection of the corre¬ 
sponding spectra show that the nuclear redshifted blob is 
associated with a strong redward asymmetry in the line 
profile, see spectrum N in Fig.l. The blueshifted region 
north-west of the nucleus persists in maps obtained with 
multi-Gaussian fits and is discussed in §5.2. 

The velocity dispersion map shows very clearly that 
the strongest line broadening occurs in the nuclear region 
and is roughly centered on the nucleus. A narrow band 
of high velocity dispersion is also observed either side 
of the nucleus, extending roughly 2" to the south-east, 
3" to the north-west and 075 above and below it. This 
band is approximately perpendicular to the axis of the 
extended line emission, which is brightest in two lobes 
that are observed north and south of the nucleus (Fig. 5). 
The greatest broadening, approximately 300 km s“^, is 
observed ~ 075 east of the nucleus. The bright lobes 
of line emission coincide with regions of relatively low 
velocity dispersion with the smallest values (a 50 km 
s“^) occurring approximately along the direction north- 
south (roughly coincident with the lobes axis). 

Gentroid velocities and velocity dispersions for the 
narrow and the broad components derived from multi- 
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Figure 4. Spectra extracted from the subregions specified in 
Fig. 6. In order to correct for Ha absorption, a stellar template 
has been fitted and subtracted from subregion D. 


Gaussian fits to the [Nil] profiles are shown in Fig. 6. The 
centroid velocity, velocity dispersion and flux maps for 
the narrow component (a 90 km s“^) are derived from 
a combination of the three-, two- and single-Gaussian 
fits. The velocity map is very similar to that derived 
from the single-Gaussian fit alone, the main difference 
being the absence of the redshifted blob at the nucleus. 
The blueshifted region north-west of the nucleus is still 
apparent. 

The centroid velocity, velocity dispersion and flux 
maps for the broad components {a ~ 200 km s“^) are 
derived entirely from the three-Gaussian fit^^. The red- 
and blue-shifted broad components are approximately 
co-spatial and dominate the line emission within the cen¬ 
tral ~ 1". However, it is notable that the brightest 
emission in the red-shifted broad component is shifted 
slightly (~ 0'.'5) north of the continuum brightness peak. 
The blue-shifted broad component is weaker overall, with 
the brightest emission slightly south-west 0'.'5) of the 
continuum peak. Evidence for line splitting (in the sense 
that two broad components are required to fit the line 
profiles) in scattered spaxels extends well beyond the 
bright central region. In fact, weak double broad com¬ 
ponents can be found out to a distance of about 2" east 
and west of the nucleus, that is, in a similar orientation 
to that of the region of enhanced velocity dispersion de¬ 
rived from the single component fit. 

The Ha, [SH] and [01] lines are characterized by the 
same features observed in [Nil], i.e. two red- and blue- 
shifted broad components concentrated in the nuclear 
region and strong narrow-lines along the north-south di¬ 
rection. 

The high velocity dispersion present in the narrow 
component map along the direction SE-NW, around the 
points (-T/5, 0'.'5) and (2'.'5, 0") in the central panel of 
Eig.6, can be attributed to the blending of the broad and 
narrow components all the way to the edge of the EOV. 
Visual inspection, see for example spectrum G in Eig.l, 
suggests a substantial decrease in the intensity of the nar¬ 
row core with respect to the broad base. Therefore, in 
that region, a single Gaussian fit would be mainly tracing 
the velocity dispersion of the broad base. Despite a num¬ 
ber of experiments, it was not possible to perform a fit 
producing a smooth transition in the velocity dispersion 
map and a perfect separation of the components. 

Velocity and velocity dispersion maps for the [EeVH 
line are presented in Eig.7. As in the case of the [Nil 
broad components, the blue- and red-shifted components 
of the [EeVH] line are approximately co-spatial, and the 
emission is concentrated within the inner 1". Erom vi¬ 
sual inspection of the datacube and the fits performed 
on individual spaxels, we consider that there is marginal 
evidence that the [EeVH] emission is elongated along the 
north-south direction. However the S/N beyond the cen¬ 
tral «0'.'5 is too low to allow a robust fit to the lines. 

4.1.2. Uncertainties 

Velocity errors for the narrow component of the [Nil] 
Ha complex are highest in regions where two component 

The broad component used in the two-Gaussian fit is not used; 
it may represent blending of red- and blue-shifted broad compo¬ 
nents. 
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Figure 5. Velocity [km s“^], velocity dispersion [km s“^] and flux [10“^® erg/s/cm^/spaxel] as derived by fitting a single Gaussian to 
the [NIIJA6583 emission line. From the flux map it is evident that the strongest emission is mainly concentrated in an unresolved region 
within 1" from the nucleus and along the north-south direction. Low-level emission is nevertheless present over the entire FOV. A systemic 
velocity Vsys = 786 km s“^ was derived from the modeling of the [Nil] kinematics and was subtracted from the velocity map. The magenta 
cross marks the continuum peak. 


fits were used: values rauge from 2 km s“^ iu the re- 
giou where the uarrow compoueut is stroug aud a broad 
base is clearly preseut, to > 50 km s~^ iu a few regious 
where the broad base is weak. Errors are siguificautly 
smaller iu regious where a siugle Gaussiau fit was used. 
Iu particular, low values (2-5 km s“^) are fouud aloug the 
directiou uorth-south, where the emissiou iu the uarrow 
compoueut is very stroug aud hues are uarrow aud siugle 
peaked. Larger values (up to 20 km s“^) are characteris¬ 
tic of the outer regious, where the S/N decreases. Errors 
iu the velocity dispersiou of the uarrow compoueut have 
a similar behavior but the rauge of values obtaiued from 
the simulatiou is smaller: typical errors for the iuuer 2" 
vary iu the rauge 10-20 km s“^. 

Errors iu the velocity of the broad compoueut iucrease 
from the ceuter outward, with the ceutral arcsecoud char¬ 
acterized by values close to 30 km s“^ iucreasiug, at the 
very edge of the fitted regiou (approximately 2" from the 
uucleus), up to 150 km s“^ for the blue compoueut aud 
close to 80 km s“^ for the red compoueut. The differeuce 
is due to the lower S/N of the blue broad compoueut. 

Errors derived for the [SII] velocity aud velocity disper¬ 
siou are very similar to those obtaiued for Ha aud [Nil]. 
Ouly two differeuces are evideut: the first is that errors 
iu the velocity of the blue broad compoueut do uot show 
the usual radial gradieut, probably because of the weak- 
uess of the hue. The secoud is that errors iu the velocity 
dispersiou of the red broad compoueut are smaller thau 
the correspoudiug errors derived from the [Nil] with val¬ 
ues raugiug betweeu 5 aud 35 km s“^ aud typical values 
of order 20 km s“^ (for comparisou, errors derived from 
the [Nil] hue rauge betweeu 5 aud 60 km s“^ showiug a 
steeper radial gradieut iucreasiug outward). 

The [01] velocity aud velocity dispersiou maps are af¬ 
fected by errors of order 1 km s“^ iu the ceuter of the 
south-westeru lobe, where the S/N is high, iucreasiug to 
10 km s“^ iu the outskirt. The uucleus, where the S/N 
for [01] is ofteu uot euough to allow a three compoueut 
fit, is characterized by errors of order 20 km s“^. The 


uorth-easteru lobe, where the sigual is lower thau iu the 
south-westeru lobe, has errors raugiug from 4 km s“^, iu 
the ceutral regiou, to 20 km s“^ aloug the rim of the lobe 
where S/N decreases. 

4.1.3. Channel maps 

Ghauuel maps were extracted iu velocity bius of 50 km 
s“^ aloug the [Nil] profile after subtractiug a systemic ve¬ 
locity Vsys = 786 km s“^ derived from kiuematical mod- 
eliug of the [Nil] velocity field. Iu order to miuimize Ha 
coutamiuatiou, uegative velocities were derived from the 
blue side of the [NII]A6548 hue, whilst positive velocities 
were derived from the red side of the [NII]A6583 hue. 
The maps are showu iu Eig.8. 

As [NII]A6548 is a factor of three weaker thau 
[NII]A6583, aud giveu that the broad blueshifted compo¬ 
ueut is iutriusically weaker thau the the redshifted oue, 
the blueshifted emissiou cauuot be mapped beyoud the 
—400 km s“^ chauuel, whereas the redshifted emissiou 
cau be mapped to +800 km s“^. 

Chauuel maps at high velocities < —200 km s“^ aud 
> 400 km s“^ are clearly domiuated by circum-uuclear 
emissiou. luterestiugly, the emissiou iu these chauuel 
maps is co-spatial despite the wide rauge of velocities 
iuvolved. Eor velocities withiu the rauge —200 : 400 km 
s~^ the hue emissiou seems to result from at least two 
kiuematic features: (i) the stroug emissiou at the uucleus 
which is isolated at extreme velocities aud is due to the 
broad spectral compoueuts, (ii) a coutributiou from the 
uarrow compoueut associated with a large-scale rotatiug 
disk. The large-scale coutributiou takes the form of off- 
ceuter lobes of emissiou which are evideut from v > —200 
km s“^ iu the south to v < 300 km s“^ iu the uorth. 

4.2. Stellar kinematics 
4.2.1. Velocity and velocity dispersion 

Iu order to derive the stellar kiuematics we fitted the 
stellar coutiuuum aud absorptiou hues usiug the Peual- 
ized Pixel Eittiug techuique (pPEX, Cappellari & Em- 
sellem 2004) aud a subset of the Vazdekis (1999) stellar 
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Figure 6. [NIIJA6583 maps. Left: velocity [km s“^]. Middle: velocity dispersion [km s“^]. Right: flux [10“^® erg/s/cm^/spaxel]. Top: 
redshifted broad component. Center: narrow component. Bottom: blueshifted broad component. Though the narrow component flux 
map shows the brightest emission to be localized along the north-south direction, low-level emission is present over the entire FOV. Given 
our fitting constraints, the velocity dispersion maps of the blue- and red-shifted broad components are the same. Labelled boxes represent 
the subregions that were rebinned to extract the representative spectra shown in Fig. 4. As evident from the bottom left panel, subregion 
N is representative of the spaxels where the broad components dominate. A systemic velocity Vsys = 786 km s“^ was derived from the 
modeling of the [Nil] kinematics and was subtracted from the velocity maps. Maps in the central row are a composite of the narrow 
component derived from three-, two- and single-Gaussian fit. Maps corresponding to the broad component of the two-Gaussians fit are not 
shown as the broad component in this case is believed to represent a blend of the blue- and red-shifted broad components included in the 
three-Gaussians fit. The broad component was included in the two-Gaussians fit only to permit accurate fitting of the narrow component. 
The black cross at (0,0) marks the continuum peak. 
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Figure 7. [FeVII]A6087 maps. Left: velocity [km s“^]. Middle: velocity dispersion [km s“^]. Right: flux [10“^® erg/s/cm^/spaxel]. Top: 
redshifted component. Bottom: blueshifted component. Note that velocity dispersions are the same by assumption. A systemic velocity 
Vsys = 786 km s“^ was derived from the modeling of the [Nil] kinematics and was subtracted from the velocity maps. 
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Figure 8. Channel maps derived from the [Nil] emission lines. Corresponding velocities are shown in the top right corner of each map 
in km s“^. Flux units are Log 10“^® erg/s/cm^/spaxel. To avoid contamination from Ha, maps corresponding to negative velocities 
have been derived from the blue side of the [NII]A6548 line, whilst maps with positive velocities have been derived from the red side of 
[NII]A6583 line. A systemic velocity Vsys = 786 km s“^ was derived from the modeling of the [Nil] kinematics and subtracted from the 
maps. Contours are equally spaced flux contours. 
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templates that were computed for an old stellar popula¬ 
tion with spectral type in the range 0-M and metallicity 
-0.7 < [Fe/H] < +0.2. Stellar centroid velocity and veloc¬ 
ity dispersions were derived, for each spaxel, after mask¬ 
ing the Na I AA5890,5896 absorption lines (which are af¬ 
fected by interstellar gas), the most prominent emission 
lines, and the O 2 A6875 sky B band. Given the spectral 
resolution of the stellar templates, 1.8 A, which is ap¬ 
proximately the same as that of our observations, 1.9 A 
at 6600 A, we did not apply any additional correction to 
the velocity dispersion. 

Although many pixels north of the nucleus were flagged 
during the fit we can see that, similarly to the gas kine¬ 
matics, the stellar velocity field (left panel in Fig. 9) dis¬ 
plays a rotation pattern in which the north-eastern side 
is receding. The stellar velocity dispersion (right panel in 
Fig. 9) reaches values close to 250 km s~^ near the center 
decreasing to 100 km s“^ toward the edges of the FOV. 
The small systematic increase, evident from south-east 
to north-west, is an artifact due to the decreasing S/N. 
The mean velocity dispersion within circles of increasing 
radii is fairly constant with a ~ 195 km s“^ within an 
aperture of radius 275 (190 pc). 

Previous estimates of the central stellar velocity dis¬ 
persion have been determined by Dressier & Sandage 
(1983), Nelson & Whittle (1995) and Garcia-Rissmann 
et al. (2005). Dressier & Sandage obtained a = 195 ± 16 
km s“^ from spectra taken with the du Pont Las Gam- 
panas 2.5 m spectrographs in the spectral region 4800- 
5400 A using a 2"x4" slit (TWHAl^^ 180 km s-^). 
This value is consistent with our findings. 

With a similar set up (aperture of 272x376 and optical 
spectra including the Mg b A5175 absorption line). Nel¬ 
son & Whittle used the Royal Greenwich Observatory 
Spectrograph on the 3.9 m Anglo Australian Telescope 
(FWHM ^ 80 km s“^) obtaining a significantly smaller 
value of cr = 120 ± 30 km s“^. Garcia-Rissmann et al. 
used the Richey-Cretien spectrograph on the 4m Mayall 
telescope at Kitt Peak National Observatory with a slit 
width of 175. From the calcium triplet they derived the 
stellar velocity dispersion using two methods obtaining 
(Ji = 123 ± 3 km s“^ and <72 = 133 ± 3 km s“^. 

The origin of the discrepancy amongst the set of values 
summarized here is unclear. 

4.2.2. Uncertainties 

Typical uncertainties in the stellar velocity field range 
from 10 km s“^, in the nuclear region, up to 15 km s“^, in 
the outer part of the field. Errors on the stellar velocity 
dispersion are slightly larger, ranging from 15 km s“^ 
near the nucleus, up to 25 km s~^ toward the outskirts 
of the FOV. 

4.3. Line fluxes 

The spatial distribution of the integrated [NIIJA6583 
line flux obtained from the single Gaussian fit to the [Nil] 
line over the whole FOV is shown in the right panel of 
Fig. 5. The strongest emission is observed within a circu¬ 
lar region of radius r ~ 1" around the nucleus. The other 
main features are two elongated regions (the “lobes”) ex¬ 
tending about 3" north and south of the nucleus in a 


direction similar to that of the compact radio emission 
reported by Nagar et al. (1999) and Mundell et al. (2009). 
Weaker [Nil] emission is present over the entire FOV at 
a level of 3 x 10“^^ erg s“^ cm“^ spaxeR^. 

The spatial distribution of the integrated line flux ob¬ 
tained from the multi-Gaussian fit is shown in the right 
column of Fig. 6 for the broad and narrow components of 
the [NII]A6583 emission line. 

As evident from the figures, the multi-Gaussian de¬ 
composition shows that the narrow component originates 
from the lobes while the broad component is essentially 
confined to the nucleus. 

Allowing for the lower spatial resolution of our data, 
the spatial distribution of the [Nil] flux is broadly con¬ 
sistent with the emission line morphologies observed in 
HST WFPG2 narrow band images in Ha + [Nil] and 
[OIII]A5007 (Ferruit et al. 2000, see also Fig. 10). The 
HST images, however, resolve substructure in the lobes, 
revealing several bright knots embedded in more diffuse 
structures, with the two brightest knots being located 
approximately 075 north and 1" south of the nucleus, re¬ 
spectively. In the HST images, the [OHI] emission has 
approximately same distribution as the Ha + [Nil] emis¬ 
sion, although the knots are more prominent in [OHI] 
whereas the more diffuse structures are more prominent 
in Ha + [Nil]. 

The HST images also reveal a “plume” of emission 
(present in both Ha + [Nil] and [OHI]) that extends 
approximately 1" to the east and north of nucleus, and 
which appears to be distinct from the north-south lobes. 
This feature has no clear counterpart in our flux maps, 
but it corresponds to a handful of spaxels characterized 
by extreme line-blending, and it may be related to the 
region of enhanced velocity dispersion seen in Fig. 5. In¬ 
deed, in the version of the Ha + [Nil] image presented by 
Ferruit et al. in their Fig. 5, a low surface brightness bar- 
like structure is present, crossing the nucleus from east to 
west, perpendicular to the lobes, and extending approx¬ 
imately 2" on each side. The plume may represent the 
inner part of this structure, which is approximately co- 
spatial with the region of enhanced velocity dispersion, 
and the velocity residuals discussed below in Section 5.2. 

The bright central region of the Ha + [Nil] flux map is 
also resolved by the HST images, into an elongated struc¬ 
ture approximately 075 long and oriented slightly west 
of north (see Fig.4 in Ferruit et ah, the feature is also 
visible in the structure map in Fig.l). Glose inspection 
shows that it is composed of two distinct components, 
to the north and south, respectively, which have higher 
excitation (as measured by the [OHI]/H(a+[NH] line ra¬ 
tio) than the surrounding more diffuse emission. These 
HST components appear to be partially resolved in the 
flux maps derived from the broad kinematic components, 
corresponding to the asymmetric brightness distribution, 
relative to the continuum peak, of the broad red- and 
blue-shifted components of the line profiles. 

The flux distribution of the blue- and red-shifted com¬ 
ponents of the [FeVH]A6087 emission line is shown in the 
right column of Fig. 7. As for the broad components of 
other lines, the redshifted component is brighter than the 
blueshifted one. 


Full width at half maximum = 2.36 cr. 


4.4. Gas exeitation 
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Figure 9. Left: velocity map as derived from the stellar absorption features excluding the [NaD] doublet. Right: stellar velocity dispersion. 
Units: km s“^. A systemic velocity r’sys,* = 743 km s“^, as derived from modeling of the stellar kinematics, was subtracted from the 
velocity field. 


Table 1 

LINE RATIOS AND ELECTRON DENSITIES 


Subr [NII]6583/Ha [OI]6366/Ha [SII](67I6 + 673I)/Ha [SII]67I6/673I Ue cm-^ 


narrow 


A 

2.45 


0.58 

0.45 


0.09 

1.07 


0.18 

1.23 ± 0.18 

200l^° 

B 

1.78 


0.01 

0.22 


0.02 

0.80 



1.24 ± 0.01 


C 

1.55 


0.06 

0.46 


0.03 

0.63 


0.2 

1.2 ± 0.04 

234t®° 

D 

1.43 


0.09 

0.36 


0.05 

0.65 


0.11 

1.6 ± 0.03 

195138 

N 

1.65 


0.01 

0.13 


0.01 

0.85 


0.01 

0.96 ± 0.02 

70411^ 








blue broad 



N 

1.12 


0.01 

0.15 


0.01 

0.24 


0.02 

2.81 ± 0.8lt 









red 

broad 



N 

1.29 



0.17 



0.57 


0.01 

0.95 ± 0.02 

7521"^ 


Note. — See the velocity map derived from the narrow component in Fig. 6 for subregion IDs. Omitted errors are smaller than 0.01. 
^Exceeds high density limit. 


When we compute the [Nil]/Ha ratio we find that the 
narrow component has values close to 1.6 approximately 
within the lobes and at the nucleus while the broad com¬ 
ponents have lower values, i.e. ~ 1.1 for the broad blue 
component and ~ 1.2 for the broad red component. 

Line ratios derived from the simultaneous fit to the 
lines in the subregion spectra are given in Table 1 and 
plotted in the Baldwin-Phillips-Terlevich (BPT) dia¬ 
grams (Baldwin et al. 1981) shown in Fig. 11. As our 
spectral coverage does not include the [OIII] and H/3 
emission lines, we estimated the [OIIIJ/H/3 ratio from 
Fig. 10 in Weaver et al. (1991). Over the whole FOV, 
the line ratios are well within the AGN photoionization 
region with the spectrum of subregion C showing ratios 
typical of LINERs and all of the others falling into the 
Seyfert region. 

This is consistent with results obtained by Bennert 
et al. (2006) who performed long slit observations of NGC 


1386, along the north-south axis of the [Nil] lobes and 
the extended [OIII] emission, finding line ratios typical 
of HII regions only beyond 6" from the nucleus. 

4.5. Electron density 

Electron densities were derived from the intensity ratio 
[SII]A6716/A6731 assuming a temperature of lO^K. The 
resulting maps are presented in Eig.l2. Eor the narrow 
component, the highest densities are found along the di¬ 
rection SE-NW, in particular in the inner 1" around the 
nucleus, reaching values « 800 cm“^. Beyond this re¬ 
gion, to the north and south of the nucleus, the density 
of the emitting gas ranges between 150 and 300 cm“^. 

High densities are also evident in the broad compo¬ 
nent maps. Median values of ~ 800 cm“^ are found for 
the blue broad component and ~ 1100 cm“^ for the red 
broad component. 

Values for the electron density derived from the spec- 
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Figure 10. Left: [OIII]A5007 image from HST observations (WFPC2-PC/F502N, proposal ID: 6419, PI: A. Wilson, Permit et al. 2000). 
Center: [OIII] contours overlaid on the [Nil] redshifted broad component flux. Right: [OIII] contours overlaid on the [Nil] narrow 
component flux. The images were aligned assuming a coincidence between the peak of the [OIII] emission and the continuum peak derived 
from our observation. 


tra extracted from the subregions identified in Fig. 6 are 
given in the last column of Table 1. It was not possible 
to derive the electron density for the blueshifted broad 
component of subregion N since the measured value of 
the [SIIJA6716/A6731 ratio (2.8) is well above the high 
density limit (approximately 1.4, Osterbrock 1989); the 
high value seems to result from a poor fit due to heavy 
blending of the lines. 

Bearing in mind the different nature and modeling of 
the data, a consistent decrease in the density as a func¬ 
tion of the distance from the nucleus was also derived 
from VLT-FORSl long slit spectroscopy by Bennert et al. 
(2006). 

5. DISCUSSION 
5.1. The stellar veloeity field 

We fitted the stellar velocity field with a kinematic 
model describing circular orbits in a plane (van der Kruit 
& Allen 1978; Bertola et al. 1991): 

Lmod{L^ Ip) = I^sysT 

Ar cos{tp — tpo) sin^ cos^^ 

{r^[sin^('0 — fio) + cos^^ cos‘^{tp — fio)] + CqCos^^}^/^ 

( 1 ) 

This yields a velocity curve that increases linearly at 
small radii and becomes proportional to at large 

radii. The parameter Vsys is the systemic velocity, A is 
the amplitude of the rotation curve, r and are the 
radial and angular coordinates of a given pixel in the 
plane of the sky, fio is the position angle of the line of 
nodes (measured with respect to the x-axis of the image 
frame, increasing counterclockwise; fio — 65° gives the 
value in degrees east of north), 0 is the disk inclination 
{0 = 0 for face-on disks). The parameter p measures the 
slope of the rotation curve where it flattens, in the outer 
region of the galaxy, varying in the range 1-1.5. Finally, 
Co is a concentration parameter, which gives the radius 
at which the velocity reaches 70% of the amplitude A. 

Because of the limited size of the region analyzed, the 


parameter p is poorly constrained by the data, therefore 
we assume p = 1, which corresponds to an asymptotically 
flat rotation curve at large radii. The stellar velocity 
curve already seems to flatten at r ~ 1" supporting this 
choice for p. The inclination 0 is also poorly constrained. 
Assuming that the gas lies on a thin circular disk, we 
adopted the value = 65° given by NED and derived 
from the 2MASS K band axis ratio of the large-scale 
structure of the galaxy. The same value is also recovered 
from PAH emission within 15" of the nucleus (Ruschel- 
Dutra et al. 2014). 

We used a Levenberg-Marquardt least-squares algo¬ 
rithm to fit the rotation model to the velocity map with 
the fit determining values for the parameters (A, Cgys, 

Co) and the center of the rotation field (xO, yO). These 
values are listed in Table 2 where results are reported for 
a fit in which the kinematical center is fixed at the con¬ 
tinuum peak and for a fit in which the center is left a free 
parameter. Maps of the observed and best fit model ve¬ 
locity field are shown in Fig. 13 together with the residual 
map for the fit performed with a fixed center. 

The fit performed with the kinematical center (xo,yo) 
left as a free parameter returns a shift from the photo¬ 
metric nucleus which is 072 north and 071 west of the con¬ 
tinuum peak. Considering our angular resolution (079), 
this is negligible and we conclude that the stellar kine¬ 
matical center coincides (within the uncertainties) with 
the photometric nucleus (N). 

The sparsity of data to the north of the nucleus in the 
velocity map (Fig. 13) might introduce a bias in the po¬ 
sition angle of the line of nodes, - 00 , toward large values 
(eastward of north). To test this possibility, we repeated 
the fit to a symmetrized version of stellar velocity map, 
which was created by reflecting the south-eastern side 
of the map about the vertical axis passing through the 
center of the image. All of the fitted parameters, and 
in particular - 00 , showed negligible variations suggesting 
that the recovered value of is not affected by the spa¬ 
tial distribution of the flagged pixels. 

5.1.1. Stability of the solution 

Dependence on the noise: to estimate the uncer¬ 
tainties in the recovered parameters due to the observed 
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Figure 11. BPT diagrams for the subregions specified in Fig. 6. 
Values for the ratio [OIII] A5007/H^ have been estimated from 
Fig.10 in Weaver et al. (1991). Vertical error bars represent the 
range of values measured by Weaver et al. within the regions of 
interest. Typical errors along the abscissa are smaller than the 
symbol size with the exception of points from subregion A and C. 
The dashed boundary lines are taken from Kewley et al. (2006). 


Table 2 

FITS TO THE OBSERVED STELLAR VELOCITY FIELD 


Parameter Notes Initial guess 


Fixed center 


A [km s-i] . 118 ± 1 100:300 

'Vsys,~k [km s“^] ... 743 ± 1 geocentric^ 500:900 

'ipo [deg] . 103 ± 0.2 60:135 

Co [arcsec] . 1.02 ± 0.03 0:5 

p . 1 fixed 

6 [deg] . 65 fixed 

xq [arcsec] . 0 fixed at N 

yo [arcsec] . 0 fixed at N 

Free center 


A [km s-i] . 117 ± 1 100:300 

'i^sys,* [km s“^] ... 744 ± 1 geocentric 500:900 

'ipo [deg] . 104 ± 0.2 60:135 

Co [arcsec] . 0.85 ± 0.04 0:5 

p . 1 fixed 

6 [deg] . 65 fixed 

xo [arcsec] . + 0.2 ^ from N -1:1 

yo [arcsec] . + 0.1 ^ from N -1:1 


Note. — ^ The formal uncertainty estimated for this coefficient 


is much smaller than 07l. ^ The heliocentric systemic velocity is 
^sys,* = 754 km s“^. 

noise level, we fitted the velocity maps recovered from the 
Monte Carlo simulation described in section §3.1. Stan¬ 
dard deviations for the distribution of the recovered pa¬ 
rameters were computed and are given in Table 2. The 
estimated uncertainties for the amplitude of the velocity 
curve, A, and the systemic velocity, r^sys, were smaller 
than 1 km s“^. 

Dependence on the initial guess: to test the de¬ 
pendence of the recovered parameters on the initial guess, 
we performed 200 fits of the stellar velocity map where 
the initial guess for each parameter was randomly drawn 
from a uniform distribution. To avoid unreasonable ini¬ 
tial guesses that would generate obviously bad fits, the 
limits of the initial guess distribution for each parameter 
were determined from visual inspection of the map to be 
fitted. The range of values allowed for the initial guess 
is given in the rightmost column of Table 2. 

The procedure was performed with both the center as 
a free parameter and the center fixed at the nucleus. In 
both cases all the fits converged exactly to the same solu¬ 
tion and none of the coefficients reached the boundaries 
allowed for the fit. 

Dependence on the value adopted for p: when 
p is fixed at values larger than 1 (implying a velocity 
curve which decreases at large radii), the amplitude in¬ 
creases steadily, reaching a value of 300 km s“^ when 
p = 1.5. A similar increase is found for the parameter 
Co (as expected from its definition) while the other fitted 
parameters, Ygys and are almost insensitive to the 
value adopted for p. 

5.2. The gas velocity field 

The velocity map derived from the narrow component 
covers the whole FOV and shows a pattern typical of 
rotation plus some distortion (see Fig. 6, first column. 
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Figure 12. Electron density maps derived from the ratio [SIIJA6716/A6731. Left: blueshifted broad component. Center: narrow 
component. Right: redshifted broad component. Units: cm“^. The color scale is arbitrarily adjusted to highlight the features of the maps. 


second row). In order to isolate non-rotational motions 
and also to determine the properties of the rotation, we 
fitted the kinematic model described by eq.l to the gas 
velocity field. The model, in this case, represents gas 
rotating in a plane with circular orbits. 

Pixels with errors in excess of 50 km s“^ were masked 
prior to performing the fit. The derived parameters are 
given in Table 3 for three realizations of the model: in 
the first case, model 1GK^ the position angle of the line 
of nodes was fixed at the value recovered from the stel¬ 
lar kinematics and the center was fixed at the continuum 
peak (which coincides with the kinematical center de¬ 
rived from the stellar velocity field); in the second case, 
model 2GK^ the position angle of the line of nodes was a 
free parameter and the center was fixed; in the third re¬ 
alization, model 3GK, both the position angle of the line 
of nodes and the center were free to vary. As for the fit 
to the stellar kinematics, for all models the parameters p 
and 0 were held fixed at values of 1 and 65° respectively. 

Assuming p = 1 implies a fiat rotation curve at large 
radii. Schulz & Henkel (2003) obtained single slit obser¬ 
vations of NGC 1386, and used the Ra and [NIIJA6583 
emission lines to derive the rotation curve out to 15" from 
the nucleus, along the galaxy major axis. From their 
Fig.4 it is evident that the rotation curve is still rising 
within the region probed in our analysis. It is therefore 
clear that neither p nor the amplitude. A, can be well 
constrained by our data. Nevertheless, the main goal of 
fitting a rotation curve to the observed velocity field is to 
isolate radial motions close to the nucleus, rather than to 
characterize the galaxy-scale rotation curve (which would 
require data from a wider range of radii). 

In model 1GK the amplitude reaches a high value (308 
km s“^) which is strongly dependent on the position an¬ 
gle assumed for the line of nodes, e.g. a variation of 1° 
in the adopted value for corresponds to a variation of 
10 - 15 km s“^ in A. 

The systemic velocity derived from this model is about 
40 km s“^ higher than the value obtained from the stel¬ 
lar kinematics. This difference could be due to the pres¬ 
ence of gaseous radial motions superposed on the rotation 


pattern. The channel maps in Fig. 8 seem to suggest a 
systemic velocity, Vcm^ which is about 70 km s“^ higher 
than derived from the modeling of gas kinematics. This 
would make the channel maps kinematically symmetric 
with respect to the zero velocity map. Nevertheless, Vcm 
is clearly not consistent with the double peaked velocity 
distribution of the narrow component (i.e. the distribu¬ 
tion of the velocity values is not symmetric with respect 
to Vqjyi). 

Models 2GK and 3GK produce relatively small resid¬ 
uals (< 30kms“^) over most of the FOV (with the ex¬ 
ceptions noted below) indicating that the gas velocity 
field is consistent with a rotating disk at a position an¬ 
gle = 88° (middle and bottom rows of Fig. 14). This 
differs by 15° with respect to the position angle obtained 
from the fit to the stellar kinematics. Offsets of compa¬ 
rable size have been found in several other Seyfert galax¬ 
ies and appear to be more common in Seyfert than in 
quiescent galaxies, perhaps indicating that weak non- 
axisymmetric perturbations are more prevalent in the 
former (Dumas et al. 2007). 

In model 3GK, where the coordinates of the kinemat¬ 
ical center are free parameters in the fit, the recovered 
center is offset 1" north of the nucleus. Offsets of similar 
size, although in different directions, were also found by 
Weaver et al. (1991), Schulz & Henkel (2003) and Ben- 
nert et al. (2006). It is possible that this offset is related 
to the distortion in the velocity field that is clearly visible 
as the blueshifted residuals (discussed below) north-west 
of the nucleus. The distortion is not due to misfits re¬ 
lated to the presence of the two broader components, it 
seems to be a genuine feature of the narrow component. 

In model 2GK, we assume that the kinematical cen¬ 
ter is that recovered from the fit to the stellar velocity 
field, which is consistent with the continuum peak, and 
hence the fit was performed with the center fixed at the 
continuum peak, (0,0). It is likely that the gas kine¬ 
matical center is affected by the presence of inflows or 
outflows, while the stellar kinematical center should be a 
more reliable tracer of the true gravitational kinematical 
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center. Therefore, we consider the parameters recovered 
from model 2GK to be more reliable than those derived 
from model 3GK. 

The overall pattern of the residuals obtained from 
models 2GK and 3GK is similar to that produced by 
model IGK. However, there are differences in the rel¬ 
ative amplitudes in the sense that for 2GK and 3GK, 
residuals in the outer regions (beyond 2" from the nu¬ 
cleus) are weaker relative to those in the inner (within 
2") regions. An interesting feature in the residual maps 
obtained from these models is a blueshifted region 2" - 
3" to the north-west of the nucleus, where the residual 
velocities reach amplitudes of 50 - 70 km s“^. This region 
is most prominent in model 2GK, but is clearly present 
in model IGK. There is evidence of a redshifted coun¬ 
terpart to the east of the nucleus. Unfortunately, the 
Ha -1- [Nil] lines are heavily blended over much of this 
region (black pixels) preventing a satisfactory line pro¬ 
file decomposition (Section 3). Nevertheless, redshifted 
residuals are present around the edges of the blended 
area in all models, strongly suggesting that a redshifted 
counterpart is indeed present. The combination of blue 
and red-shifted residuals together indicate the present of 
a distinct kinematical component spanning the nucleus, 
approximately along the galaxy minor axis. 

Comparison with the velocity dispersion map pre¬ 
sented in Fig. 5 indicates that this residual is associ¬ 
ated with the region of high velocity dispersion extend¬ 
ing from —2" to +3" along the SE-NW direction (i.e., 
roughly along the galaxy minor axis) and smoothly tran¬ 
sitioning to lower velocity dispersions. As noted in Sec¬ 
tion 4.3, these kinematic features are approximately co- 
spatial with a low surface brightness bar-like structure 
seen in the HST Ha + [Nil] image of Ferruit et al. 2000. 

The velocity dispersion map presented in Fig. 5 shows 
that these residuals are associated with a region of high 
velocity dispersion extending from —2" to +3" along the 
SE-NW direction and smoothly transitioning to lower 
velocity dispersion values. 

As evident from the velocity maps in Fig. 6 the two 
broad components are systematically blue- and red¬ 
shifted with respect to the larger-scale narrow compo¬ 
nent, suggesting the presence of either an inflow or an 
outflow. Their mean velocities relative to the systemic 
velocity derived from the kinematics of the narrow com¬ 
ponent are —140 km s“^ and 250 km s“^ respectively. 
These broad components appear co-spatial in our flux 
maps, but it seems reasonable to identify at least the 
central part of the region emitting them with the two 
components (northern and southern) of the bright cen¬ 
tral feature in the HST emission line images of Ferruit 
et al. (2000). However, evidence for line splitting, in the 
sense that two broad components are required to fit the 
line profiles, extends well beyond the bright, unresolved 
central region. In fact, weak double broad components 
can be found out to a distance of approximately 2" east 
and west of the nucleus, that is, in a similar orientation 
to that of the region of enhanced velocity dispersion de¬ 
rived from the single-component fit. 

The [FeVH] line profile is also dominated by blue- 
and red-shifted components, although with more extreme 
mean velocities (—250 and 270 km s“^. Fig. 7). This line 
originates closer to the AGN, therefore velocities higher 


than those derived for the low-ionization species are ex¬ 
pected. The emission in this line is also centrally con¬ 
centrated (within 1" of the nucleus) and again, it seems 
reasonable to identify the blue- and red-shifted kinematic 
components with the central structures visible in the 
HST image. 

5.2.1. Stability of the solution 

Dependence on the noise: as previously described 
for the stellar kinematics, we fitted the velocity maps 
recovered from the Monte Carlo simulation described in 
section §3.1. Standard deviations for the distribution of 
the recovered parameters were computed and are given 
in Table 3 as an estimate of the uncertainty. 

Dependence on the initial guess: as for the stel¬ 
lar kinematics, we tested the dependence of the recov¬ 
ered parameters on the initial guesses and again we find 
that using initial guesses drawn randomly from the range 
specified in Table 3 the solution is always stable converg¬ 
ing to the same values. 

Dependence on the value adopted for p: the am¬ 
plitude, A, and the concentration parameter, cq, are the 
most sensitive to the value adopted for p. We assumed 
p = 1, which is justified by the flattening of the rotation 
curve at large radii (Schulz & Henkel 2003). When p is 
larger than 1, the amplitude immediately reaches high 
values (i.e. larger than 400 km s“^ for model IGK and 
larger than 280 km s“^ for models IGK and 2GK). The 
systemic velocity and are almost insensitive to p. 

5.2.2. Comparison with previous work 

Although parameters like the amplitude. A, and cq 
may not be reliable due to the limited spatial extent 
of our data, we believe that the other parameters are 
reliably and accurately determined. Perhaps the most 
important of these parameters is the systemic velocity. 
Recall that the systemic velocity derived from stellar 
kinematics should be the most reliable measure of the 
galaxy recessional velocity as the value derived from gas 
kinematics may be affected by non-rotational motions, 
such as inflows or outflows. 

Various estimates of the systemic velocity have been 
made at lower spatial resolution using the gas kine¬ 
matics derived from long-slit observations of the Ha 
and [NH]A6583 emission lines (e.g. Weaver et al. 1991; 
Storchi-Bergmann et al. 1996; Schulz & Henkel 2003). 
All of these measurements produce values in the range 
868-890 km s“^, while the heliocentric systemic velocity 
derived here is Vsys gas = 797 ± 1 km s“^ (from model 
2GK). 

A contribution to the discrepancy comes from our de¬ 
cision to fix the gaseous kinematical center at the stellar 
kinematical center. Indeed, the discrepancy decreases 
when the center is fitted as a free parameter, but there 
is still a difference Av > 50 km s“^ with previous mea¬ 
surements. In addition, the lower resolution of previous 
studies would have made it more difficult to disentangle 
the rotating velocity field from the additional compo¬ 
nents. However it should be noted that Schulz & Henkel 
derive a systemic velocity of 877 km s“^ relying solely on 
the outer portion of their velocity curve, i.e. using data 
between S" and 15", where the curve flattens. When 
measuring the velocity at the brightness peak, they ob- 
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Table 3 

FITS TO THE OBSERVED GASEOUS VELOCITY FIELD 


Parameter Notes Initial guess 


Geometry from stellar kinematics (model IGK) 


A [km s-i] . 308 ± I 0:900 

'i^sys,gas [km s“^] .. 787 ± I geocentric^ 500:900 

'ipo [deg] . 103 fixed 

Co [arcsec] . 4.93 ± 0.01 0:5 

p . 1 fixed 

6 [deg] . 65 fixed 

xq [arcsec] . 0 fixed at N 

yo [arcsec] . 0 fixed at N 

Fixed center (model 2GK) 


A [km s-i] . 226 ± 1 180:300 

Vsys,gas [km s“^] .. 786 ± 1 geocentric^ 500:900 

'ipo [deg] . 88 ± 0.3 60:100° 

Co [arcsec] . 3.08 ± 0.02 0:5 

p . 1 fixed 

6 [deg] . 65 fixed 

xo [arcsec] . 0 fixed at N 

yo [arcsec] . 0 fixed at N 

Free center (model 3GK) 


A [km s-i] . 213 ± 1 180:300 

'i^sys,gas [km s“^] .. 810 ± 1 gcocentric^ 500:900 

ipQ [deg] . 88 ± 0.1 60:100° 

CO [arcsec] . 2.57 ± 0.02 0:5 

p . 1 fixed 

6 [deg] . 65 fixed 

Xo [arcsec] . + 0.4 from N -1:1 

yo [arcsec] . + 0.9 ^ from N -1:1 


Note. — ^ The estimated uncertainty for this coefficient is much 
smaller than 07l. ^The heliocentric systemic velocity is derived by 
adding 11 km s“^ to the geocentric value. 

tain a value (t;sys = 789 ±15 km s“^) consistent with 
our findings. 


The kinematics of the coronal line [FeVII]A6087 
in NGC 1386 have previously been investigated by 
Rodriguez-Ardila et al. (2006). They performed single 
slit observations of the galaxy, with the slit aligned along 
the north-south direction. As in this work, they detected 
a double peaked emission line with the redshifted compo¬ 
nent stronger than the blueshifted one. They observed 
emission up to about 106 pc (l'/4) north and south of 
the nucleus (the spatial scale has been adjusted to our 
adopted distance). The maps presented in Fig. 7 show a 
smaller spatial extension. Visual inspection of the dat- 
acube suggests that emission is present over a scale com¬ 
parable with that proposed by Rodriguez-Ardila et ah, 
however the S/N is not high enough to allow a robust fit 
to the spectra. As proposed by Rodriguez-Ardila et ah, 
the presence of this emission, roughly aligned with the 
bright lobes visible in the flux maps, and probably asso¬ 
ciated with the bright nuclear knots revealed by Ferruit 
et al. (2000), is suggestive of a bipolar outflow along the 
radiation cone. 

6. INTERPRETATION 


Our analysis of the IFS data cube has revealed mul¬ 
tiple components of line-emitting gas characterized by 
distinct kinematic properties. In particular, we iden¬ 
tify three main kinematic components: (1) a low velocity 
dispersion (“narrow”) component that extends over the 
entire FOV, (2) a higher velocity dispersion (“broad”) 
component observed within approximately 2" of the nu¬ 
cleus, but mostly concentrated within 1" and apparently 
coincident with the two nuclear knots resolved in the 
HST observation of Ferruit et al. (2000), and (3) a re¬ 
gion of enhanced velocity dispersion, associated with red- 
and blue-shifted velocity residuals and with a faint emis¬ 
sion structure seen in the Ferruit et al. HST imaging, 
points to the presence of a distinct kinematic component 
roughly perpendicular to the galaxy disk. 

The velocity field of the narrow component is consis¬ 
tent with gas predominantly rotating in circular orbits 
with inclination consistent with the galaxy disk. How¬ 
ever, fits to the velocity field reveal residuals indicating 
that significant deviations from rotation about the axis 
of the large-scale disk are also present. 

6.1. The narrow component - disk rotation 

The presence of multiple components of emission lines 
is evident from the total flux map presented in Fig. 5, as 
well as from the flux distributions in the third column of 
Fig. 6, and the channel maps of Fig. 8. The bright lobes 
extending along the direction north-south, and the low- 
level emission present over the entire FOV, are associ¬ 
ated with the narrow kinematic component. They show 
rotational motion and low velocity dispersion. There¬ 
fore, although the double-lobed morphology of the nar¬ 
row component is suggestive of an outflow, the velocity 
field does not support this interpretation. The morphol¬ 
ogy and kinematics can be reconciled if we interpret the 
bright lobes as due to emission from photoionized gas in 
the regions of the galaxy disk that are intersected by the 
biconical AGN radiation field that is predicted by the 
AGN unified model (e.g. Antonucci 1993). Our interpre¬ 
tation is based on the model proposed for NGG 4151 by 
Robinson et al. (1994): if the edge of the radiation cone 
emerges from the disk at a shallow angle, then the in¬ 
tersection with the disk surface traces an ellipse, which 
would result in an elongated surface brightness distri¬ 
bution similar to what is observed (see Section 6.5 and 
Fig. 16). 

Fig. 10 shows that the [OHIJA5007 emission matches 
very well the flux distribution observed for the [Nil] emis¬ 
sion and the other low-ionization species included in our 
observation. It therefore seems reasonable to assume 
that the extended [OHI] emission has the same origin 
as the low-ionization lines. 

6.2. The broad component - nuclear outflow 

The second kinematic component - the broad compo¬ 
nent - is associated with the strong nuclear line emission 
responsible for the central “blob” in the flux distribu¬ 
tion maps. It is mostly confined within 1" from the nu¬ 
cleus, but shows some extension out to approximately 2". 
This component is actually double, and appears in the 
emission-line profiles as two components, one blueshifted 
and the other redshifted relative to the systemic velocity. 
They are co-spatial and have a relatively large velocity 
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Figure 13. Left: velocity map [km s as derived from the stellar absorption features excluding the [NaD] doublet. Center: rotation 
model fit. Right: residual. 


dispersion (a 220 km s“^). A reasonable interpreta¬ 
tion is that they represent the approaching and receding 
sides of an outflow. 

Spectra presented in Fig. 2 and the flux maps in the 
right column of Fig. 6 show that the red broad compo¬ 
nent is brighter than the blue. This type of asymmetry 
(redshifted gas brighter than blueshifted gas) is not what 
would usually be expected in an outflow. Indeed, emis¬ 
sion line blueshifts or blue asymmetries are usually con¬ 
sidered evidence for outflows (e.g. Zamanov et al. 2002). 
Here we suggest that the origin of the asymmetry in the 
brightness of the two components is due to the combi¬ 
nation of the dust distribution and morphology of the 
outflow: we propose that the broad components orig¬ 
inate from outflowing discrete clouds containing dust. 
These clouds are photoionized by radiation from the 
AGN around which they are distributed, which implies 
that the ionized face of each cloud is always directed 
toward the AGN. Therefore blueshifted clouds preferen¬ 
tially present to the observer their non-ionized side; on 
the other hand, the redshifted clouds present their ion¬ 
ized faces to the observer, thus appearing brighter. 

The picture outlined above also explains the asymme¬ 
try in the velocity of the broad components (flbiue ~ 
— 140 km s“^, Vred ~ 250 km s“^): blueshifted clouds 
with maximal line-of-sight velocity are those at small pro¬ 
jected distances from the AGN, e.g. cloud 1 in Fig. 15. 
According to our proposed scenario, radiation from these 
clouds is strongly attenuated as they show the observer 
their non-ionized face. Blueshifted clouds located at 
larger projected distances, and therefore characterized 
by smaller line-of-sight velocities, contribute a greater 
fraction of the blueshifted emission as their ionized face 
is partially visible, e.g. cloud 2 in Fig. 15. On the other 
side, redshifted clouds at small projected distances from 
the AGN (and therefore characterized by maximal line- 
of-sight velocities) always show the observer their ionized 
faces. As a result, the average velocity measured for the 
blueshifted clouds is smaller than the velocity measured 
for the redshifted counterparts. 

It seems reasonable to identify the blue- and red-shifted 
broad components of the line profile with the two com¬ 


ponents (northern and southern) of the bright central 
feature in the HST emission line images of Ferruit et al. 
(2000). If this association is correct, it suggests that the 
kinematic components are the approaching and receding 
sides of a bipolar outflow that is aligned with the radia¬ 
tion cone axis. 

Maps for the [FeVII] emission lines are consistent with 
those derived for the [Nil] broad component, however 
there is marginal evidence that the emission is extended 
along the radiation cones. This is consistent with the idea 
that the blue- and red-shifted broad components come 
from a bipolar outflow aligned with the AGN radiation 
cones. 


6.2.1. Mass outflow rate 

We estimate the mass outflow rate as the ratio between 
the mass of the outflowing gas and the dynamical time 
at the nucleus. The gas mass associated with the broad 
component is estimated as: 

= ( 2 ) 

where m^ is the proton mass, rig ^ 1200 cm“^ and V 
are the electron density and the volume of the region 
where broad components are detected, / is the Ailing 
factor. The Ailing factor and volume are eliminated by 
combining eq.2 with the following expression for the Ha 
luminosity: 


LH^--nlmAT)Vf (3) 

with jH„(T) = 3.534 X cm^ erg s”! at T=10000 K 

(Osterbrock 1989). Therefore the gas mass is expressed 
as: 


Mg = 


(T) ■ 


( 4 ) 


In the previous section we argued that emission from 
the blueshifted broad component is attenuated by dust 
embedded within the emitting clouds, whilst clouds re¬ 
sponsible for the redshifted broad component show the 
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Figure 14. Top: model IGK, the center is fixed at the continuum peak and the line of nodes is fixed at the value derived from stellar 
kinematics. Middle: model 2GK, the line of nodes is a free parameter and the center is fixed at the continuum peak. Bottom: model 
3GK, line of nodes and center are free parameters. Left: velocity map [km s“^] as derived from the narrow component of the [NIIJA6583 
emission line; pixels with errors larger than 50 km s“^ are masked (in black). Genter: rotation model fit [km s“^]. Right: residual [km 
s“^]. Magenta cross: kinematical center. Black cross: continuum peak. Maps are shown after subtraction of the systemic velocity derived 
from each fit. The fit parameters are given in Table 3. 
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observer their ionized side. Therefore, we assume that 
the red broad component gives a better representation of 
the intrinsic emission produced by the outflowing gas and 
we estimate the total Ha luminosity due to the outflow 
as twice the Ha luminosity from the red broad compo¬ 
nent, yielding = 2.5 x 10^^ erg s“^. The mass of the 
ionized gas is, therefore. Mg ^ 5 x 10^ Mq. This repre¬ 
sents the amount of gas responsible for the Ha emission 
from the two “broad” components. 

We estimate the dynamical time as the ratio of the ra¬ 
dius of the region where broad components are observed 
2" or 152 pc) to the mean velocity of the outflow 
(250 km s“^). This gives ^ Q x 10^ yr. The corre¬ 
sponding mass outflow rate is M ^ 0.1 Mq yr“^ (note 
that this is a lower limit as it represents the outflowing 
mass associated with only the ionized side of the clouds). 
For comparison, the mass accretion rate implied by the 
bolometric luminosity is Mace = Lboi/vc'^ - 3 X lO-^M© 
yr“^, where we assume a mass-to-energy conversion efh- 
ciency r] = 0.1, Lijoi ~ 90I/[o///] (Dumas et al. 2007) and 
L[oiii] ~ 2 X 10^^ erg s“^ from the flux measurement in 
Storchi Bergmann & Pastoriza (1989). 

6.3. Equatorial rotation or expansion? 

In general, the narrow component velocity held is dom¬ 
inated by circular rotation consistent with the rotation 
of the large-scale galactic disk. However, there is strong 
evidence for a distinct kinematic component extending 
2" - 3" either side of the nucleus, and approximately ori¬ 
ented along the minor axis of the large scale disk. This 
component appears to be associated with a low surface 
brightness bar-like structure that is seen crossing the nu¬ 
cleus east to west in the HST WFPC2 image obtained 
by Permit et al. (2000), Sec. 4.3. It is characterized by 
increased velocity dispersion and blue- and red-shifted 
velocity residuals to the west and east of the nucleus, re¬ 
spectively. Notably, the orientation of this structure is 
roughly perpendicular to the axis of the AGN radiation 
bicone. 

The velocity residuals can be interpreted as either ra¬ 
dial outflow or rotation. In the first case, outflow is fa¬ 
vored over inflow on the basis of geometrical arguments: 
the blue residual is associated with the near side of the 
galaxy disk. In the second case, the rotation is about an 
axis approximately coincident with the axis of the AGN 
radiation cone. In either case, the outflow/rotation ap¬ 
pears to be occurring in a plane inclined at large angle 
with respect to the large-scale disk and roughly perpen¬ 
dicular to the AGN radiation bicone. We propose that 
this is the equatorial plane of the torus. It is possible 
that the residuals result from a combination of rotation 
and outflow. Several models (e.g. Krolik & Begelman 
1986; Elitzur & Shlosman 2006; Dorodnitsyn et al. 2008; 
Keating et al. 2012; Wada 2012; Dorodnitsyn & Kail- 
man 2012) propose that the torus itself is formed in an 
approximately equatorial magneto-hydrodynamical wind 
or that a wind is blown off the torus. The observed resid¬ 
uals, the enhanced velocity dispersion, and the circum- 
nuclear surface brightness structure, could conceivably 
originate from an extension to larger scales (or from a 
remnant) of such a wind. A cartoon illustrating the pro¬ 
posed scenario is shown in Fig. 16. 


Vt2 V2 



Figure 15. Outflowing clouds (circle 1 and 2) are ionized by the 
AGN accretion disk (ellipse on the left). Blueshifted clouds pref¬ 
erentially show to the observer their non ionized face (in black). 
Their line emission is attenuated by dust embedded in the cloud. 
Clouds at large projected distances from the nucleus (large values 
of a, corresponding to small line-of-sight radial velocities, Vr) are 
brighter as their ionized face is partially visible to the observer. As 
a result, the average velocity measured for the blueshifted clouds 
is skewed toward small values. 

6.4. Velocity residuals in the large-scale disk 

In the outer regions of the FOV, at distances of 2" 
- 5" from the nucleus, we And blueshifted residuals in 
the E-NE arc, while redshifted residuals are found to 
the N-NW and W-SW. This system is most prominent 
in the residual maps obtained from model IGK, where 
the line of nodes and the kinematical center are fixed at 
the values determined from stellar kinematics, but can 
also be discerned at much lower amplitude in the resid¬ 
ual maps obtained from the other two models (2GK and 
3GK). The blueshifted velocities appear in the far side of 
the galaxy (the south-eastern side, see Eig.I) while the 
redshifted velocities appear in the near side, and are spa¬ 
tially coincident with the arc-like dust lanes highlighted 
in the structure map presented in Eig.I. The dust lanes 
appear to trace part of a nuclear spiral, and we suggest 
that these velocity residuals may trace gas streaming in¬ 
wards along the spiral. 

6.5. The geometry and kinematics of the nuclear region 

Here we attempt to place the various components of 
the circum-nuclear ionized gas in NGC 1386, as charac¬ 
terized by distinct kinematics and spatial distributions, 
in the context of an overall structure for the nuclear re¬ 
gion. We also consider the relationship of these com¬ 
ponents to the classical emission line regions of AGN, 
namely the narrow line region (NLR) and the extended 
narrow line region (ENLR). Our proposed configuration 
for the nuclear region of NGC 1386 is shown schemati¬ 
cally in Eig.I6. 

The gas emitting the narrow kinematic component re¬ 
sides in a disk, which rotates approximately in the plane 
of the galaxy. The line emission from this component ex¬ 
tends over the entire EOV (690 x 530 pc), but is brightest 
in a double-lobed structure that results from photoion¬ 
ization of the disk gas where the AGN radiation cones 
intersect the disk (Eig.I6, left and center). The geome¬ 
try of the radiation cones can be constrained by requiring 
(I) that the disk surface area bounded by the intersec¬ 
tion with the cones approximately matches the extension 
and orientation of the bright lobes, and (2) that the line 
of sight to the nucleus does not pass within the radia¬ 
tion cone (i.e., the line of sight must intercept the torus, 
since NGC 1386 is a Seyfert 2). With these constraints. 
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Figure 16. The proposed configuration for the nuclear region of NGC 1386 (see Section 6.5). Left: the narrow component of the emission 
lines originates from ionized gas rotating (long-dashed arrow) and infiowing (short-dashed blue and red arrows) in the plane of the galaxy 
(large grey disk). The yellow cones schematically represent the AGN radiation field which escapes the torus in bicones aligned with the 
torus axis. The radiation cone axis is inclined at 38° to the disk plane and the cones have a half-opening angle of 34°. Center: the bright 
lobes visible in the fiux maps are due to photoionization of the disk gas where the radiation cones (omitted for clarity) intersect the disk. 
Right: within the central 1", the broader emission line components arise mainly from dusty clouds (green blobs) participating in a 
bipolar radial outfiow (blue and red dashed arrows) aligned with the radiation cone axis. Blue- and red-shifted velocity residuals, which are 
associated with increased velocity dispersion, suggest that rotation (solid green arrow) and/or outfiow (short black arrows) is occurring on 
larger scales (2 - 3'') in a plane approximately perpendicular to the radiation cone axis (thin green disk). The large-scale disk associated 
with the narrow emission line component is omitted for clarity. Cartoon created using the Shape program (Steffen et al. 2010). 


we find that the observed morphology is approximately 
reproduced when the radiation cone axis is inclined at 
38° to the disk plane and the cones have a half-opening 
angle of 34° (however, these values should be considered 
illustrative rather than tightly constrained). 

The extended emission-line lobes, therefore, do not 
represent a physical (matter bounded) structure, such as 
an outflow or, as suggested by Schulz & Henkel (2003), 
an edge-on disk; they are merely those parts of a disk 
rotating in the galactic plane that are illuminated (and 
photoionized) by the AGN radiation emerging from the 
torus opening. In this respect, they are similar to the ^ 
kpc scale “extended narrow line regions” that are present 
in many Seyfert galaxies (i.e., interstellar gas clouds pho¬ 
toionized by the anisotropic AGN radiation field, Unger 
et al. 1987). The overall similarity between the mor¬ 
phologies of the extended (> 1" north and south of the 
nucleus) [OIIIJA5007 and Ra + [Nil] emission in HST im¬ 
ages (Ferruit et al. 2000) and the flux distributions in the 
narrow component of the [Nil] (and Ha) emission lines, 
derived from our IFS data, implies that the extended 
[OHI] emission has the same origin. The sub-structure 
in the “lobes” revealed by the HST images presumably 
reflects the non-uniform distribution of gas clouds in the 
disk. 

It is worth noting that Ghandra observations of NGC 
1386 presented by Li & Wang (2013) show the presence of 
elongated X-ray emission, extending over few kpc, with 
an orientation consistent with that of the radiation cones 
inferred here. 

As noted in Section 6.4, the pattern of velocity resid¬ 
uals suggests that inward spiral streaming motions may 
also be present in the disk, although the amplitudes are 
dependent on whether or not the gas disk is assumed to 
be co-planar with the stellar rotation. 


We consider the gas emitting the broad double compo¬ 
nent to be participating in a nuclear outflow, the blue- 
and red-shifted components representing the approach¬ 
ing and receding sides of the flow, respectively. Our data 
do not have sufficient spatial resolution to strongly con¬ 
strain the geometry of the outflow. However, the bright¬ 
est emission is concentrated within the central 1" and 
may be associated with the elongated double structure 
which is visible in both the structure map presented in 
Fig.l, and in the HST WFPC2 [OHI] and Ha + [Nil] im¬ 
ages of Ferruit et al. (2000). If this is the case, the mor¬ 
phology suggests a bipolar outflow approximately aligned 
with the radiation cone axis. The southern HST com¬ 
ponent would then be identified with the approaching, 
blueshifted, side and the northern component with the 
receding, redshifted side. 

The characteristic properties of the broad component 
gas, in particular its proximity to the nucleus, electron 
density (^1000 cm“^), velocity dispersion (?^200 km s“^) 
and its complex kinematics, lead us to identify it as the 
NLR in NGC 1386. In our illustrative cartoon, this 
is represented by the green blobs in the right panel of 
Fig. 16, each one representing a distribution of outflow¬ 
ing clouds. 

In addition to the two main kinematic components, 
that is, rotation in the plane of the galaxy and a nuclear 
outflow, there is strong evidence for another component 
that involves rotation and/or outflow in the equatorial 
plane of the torus, and extends to approximately 2 - 3" 
(^200 pc) either side of the nucleus. We speculate that 
this structure (represented by the green disk in the right 
panel of Fig. 16) may be evidence for the accretion-disk 
winds that have been proposed to explain the formation 
of the torus. The excitation mechanism for the relatively 
faint line emission associated with this kinematic com¬ 
ponent is unclear. The gas is (of necessity) located out- 
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side the AGN radiation cones. However, the line ratios 
in subregions C and D, which sample the emission west 
and east, respectively, of the nucleus, are generally con¬ 
sistent with AGN photoionization, although subregion G 
in particular, exhibits lower excitation characteristic of 
LINERs. Thus, it is possible that this gas is photoionized 
by an attenuated AGN radiation field “leaking” through 
a clumpy torus. Alternatively, the enhanced velocity dis¬ 
persion in this region suggests that shock ionization could 
also play a role. 

7. SUMMARY AND CONCLUSIONS 

Using the GMOS integral field unit we observed the 
inner 690 pc of the nearby Seyfert 2 spiral galaxy NGG 
1386 in the spectral range 5600 - 7000 A which in¬ 
cludes a number of prominent emission lines ([OIJA6302, 
[FeVII]A6087, [NII]AA6548,6583, Ha, [SH]AA6716,6731). 
We modeled the profiles of these lines to produce ve¬ 
locity, velocity dispersion and flux maps with a spatial 
resolution of 68 pc and a spectral resolution of 66 km 
s“^. 

The emission-line flux maps, as well as flux distribu¬ 
tions in channel maps, are characterized by two elongated 
structures to the north and south of the nucleus, a bright 
nuclear component with FWHM « O'.'Q, and low-level 
emission extending over the whole FOV. We argue that 
the elongated structures are due to emission from gas in 
the galaxy disk photoionized by AGN radiation via a ge¬ 
ometry in which the AGN radiation bicone intercepts the 
galactic disk at a shallow angle. Fainter emission consis¬ 
tent with AGN photoionization is also observed beyond 
the elongated structures suggesting that some radiation 
from the AGN escapes through the obscuring torus. 

The gas has three distinct kinematic components. 

1. A low velocity dispersion component (a ~ 90 km 
s“^), which is consistent with gas rotating in the 
galaxy disk and includes the bright elongated struc¬ 
tures (the “lobes”) observed in the flux map. There 
is no kinematic evidence that these structures are 
distinct from the galaxy disk gas; they are partici¬ 
pating in the general rotation and, as noted above, 
can be explained by partial illumination of the disk 
by the AGN radiation cones. However, subtraction 
of rotation models reveals velocity residuals in the 
outer regions of the FOV that suggest inflow, per¬ 
haps in the form of streaming along dusty nuclear 
spirals. We estimate the electron density from the 
[SH] doublet, finding that the narrow component 
has Tie ~ 1000 cm“^ at the nucleus and 200 
cm“^ in the bright extended lobes. 

2. A higher velocity dispersion component (a ~ 220 
km s“^) characterized by blueshifted and red- 
shifted sub-components at average velocities of 
— 140 and +250 km s“^, respectively. This com¬ 
ponent corresponds to the bright nuclear emission 
observed in our flux maps and is attributed to a 
compact (< 1") bipolar outflow aligned with the 
radiation cone axis. It appears to be resolved into 
two bright knots by HST observations, which can 
be identified with the approaching and receding 
sides of the outflow. We identify it as the nar¬ 
row line region in NGG 1386. The broad compo¬ 


nents have electron densities of ~ 800 cm“^ for 
the blueshifted and 1100 cm“^ for the redshifted 
components, respectively. Using the Ha luminosity 
we estimate a lower limit for the mass outflow rate 
of Mout ^0.1 Mq yr which is much larger than 
the accretion rate necessary to sustain the AGN 
bolometric luminosity, that is Macc ^ 3 x 10-4 M© 
yr-i. 

3. There is strong evidence for a third kinematical 
component in the form of velocity residuals that 
are spatially coincident with a region of enhanced 
velocity dispersion and a faint, bar-like emission 
structure revealed by HST imaging. This compo¬ 
nent extends 2"-3" either side of the nucleus, ap¬ 
proximately perpendicular to the major axis of the 
galaxy disk. The presence of blueshifted residuals 
on the near (west) side and redshifted residuals on 
the far (east) side of the galaxy is consistent with 
outflow and/or rotation in a plane that is approx¬ 
imately perpendicular to the AGN radiation cone 
axis. We speculate that this is a wind which is both 
outflowing in a plane roughly coincident with the 
equatorial plane of the torus, and rotating about 
the axis of the radiation cones. 

In some respects, the structure and kinematics of the 
circum-nuclear region of NGG 1386 are consistent with 
the “classical” AGN unification model in which a hy¬ 
drostatic, optically thick torus prevents isotropic escape 
of the UV-optical radiation from the AGN. In particu¬ 
lar, the bright extended lobes of line emission can be ex¬ 
plained by AGN photoionization where biconical beams 
of ionizing radiation intersect the galaxy disk. There is 
also good evidence for a compact bipolar outflow along 
the axis of the AGN radiation cones; such outflows ap¬ 
pear to be a common feature of Seyfert galaxies (e.g. 
Fischer et al. 2013). 

However, we also find features that are not easily 
explained in the context of the “classical” unification 
model. The emission line ratios are generally consistent 
with AGN photoionization over the whole FOV, even 
for the faint emission outside the inferred boundaries for 
the AGN radiation cones. This suggests that attenuated 
AGN ionizing radiation is able to escape at large angles 
from the radiation cone axis, favoring models in which 
the torus consists of an ensemble of discrete, optically 
thick clumps (e.g. Nenkova et al. 2008). 

Perhaps more intriguing, however, is the evidence that 
rotation and/or outflow is also occurring in the equa¬ 
torial plane of the torus. In recent years, our group 
has found evidence for the presence of similar compact 
(< 100 pc) equatorial outflows, i.e. perpendicular to 
the axis of the radio jet or to the axis of the ionization 
cones, in several other active galaxies (Gouto et al. 2013; 
Schnorr-Miiller et al. 2014b; Riffel et al. 2014); an equa¬ 
torial outflow has also been suggested by Su et al. (1996) 
to explain the features observed in the radio emission 
of NGG 5929. Perhaps we are finding evidence of torus 
winds similar to those predicted by recent models (e.g. 
Keating et al. 2012; Dorodnitsyn & Kallman 2012, and 
references therein). These observations hint that equato¬ 
rial outflows on ^ 100 pc scales may be a common feature 
of AGNs. However, to investigate their nature and ubiq- 
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uity among the AGN population, high spatial resolution 
spectroscopy of a larger sample and detailed modeling to 
determine the kinematic and physical properties of the 
gas will be required. 
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